Structure and stability of Mott-insulator shells of bosons trapped in an optical lattice 
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We consider the feasibility of creating a phase of neutral bosonic atoms in which multiple Mott- 
insulating states coexist in a shell structure and propose an experiment to spatially resolve such a 
structure. This spatially-inhomogeneous phase of bosons, arising from the interplay between the 
confining potential and the short-ranged repulsion, has been previously predicted. While the Mott- 
insulator phase has been observed in an atomic gas, the spatial structure of this phase in the presence 
of an inhomogeneous potential has not yet been directly probed. In this paper, we give a simple 
recipe for creating a structure with any desired number of shells, and explore the stability of the 
structure under typical experimental conditions. The stability analysis gives some constraints on how 
successfully these states can be employed for quantum information experiments. The experimental 
probe we propose for observing this phase exploits transitions between two species of bosons, induced 
by applying a frequency- swept, oscillatory magnetic field. We present the expected experimental 
signatures of this probe, and show that they reflect the underlying Mott configuration for large 
lattice potential depth. 

PACS numbers: 32.80.Pj,03.75.Lm 



I. INTRODUCTION 

Recent experimental results involving quantum degen- 
erate atom gases trapped in optical lattices have stimu- 
lated interest from the perspectives of condensed matter 
and quantum information science. Ultra-cold atoms con- 
fined in an optical lattice are predicted to display a rich 
variety of quantum phases (see 0,0, for example). The 
ability to precisely control the physical parameters of this 
system enables probing a vast range of related funda- 
mental physics including quantum phase transitions, the 
behavior of collective excitations in the quantum regime, 
and the physics of defects. The properties of these phases 
may permit new applications, such as neutral atom quan- 
tum computing or quantum simulation 

Neutral bosons trapped in optical lattices can exhibit 
complex spatial configurations of co-existing superfiuid 
and Mott-insulator phases 0. There have been detailed 
studies of the superfiuid to Mott-insulator phase transi- 
tion in the atomic system using techniques such as quan- 
tum Monte Carlo simulation and mean-field theory (see 
^lEIEQ for examples). In this paper we offer a practi- 
cal and simple recipe for realizing specific Mott-insulator 
configurations and we study the limitations posed by re- 
alistic experimental conditions. Motivated by proposals 
for using the Mott-insulator state to initialize a fiducial 
state for neutral atom quantum computing 8,'9'|, wc fo- 
cus on the tight-binding regime where intcr-sitc tunneling 
is suppressed. We formulate a straightforward method, 
primarily using counting arguments, for obtaining the 
constraints on lattice parameters that yield states of def- 
inite spatially varying occupation number. 

The effects of finite temperature, tunneling, and dissi- 
pation in a peri odic potential have been studied in the lit- 
erature [Tolllll IT3 | . Atomic systems are a unique tool for 
studying outstanding questions regarding these phenom- 



ena. Furthermore, deviations from the zero-temperature 
and tightly-bound ground state will have an important 
impact on potential applications using atoms trapped in 
a lattice. We investigate processes which disturb the sta- 
bility of the Mott-insulator ground state and estimate 
their influence for characteristic experimental conditions. 
Specifically, we initiate an analysis of the emergence of 
superfiuid order when there is a finite but small degree of 
tunneling and we study thermally-activated defects and 
their dynamics. Our work is not meant to be an exhaus- 
tive study of stability. Instead, we take the first steps 
toward addressing these considerations and point to ex- 
isting literature from the condensed matter community 
appropriate to the atomic system. 

To complement our theoretical work, we propose a set 
of experiments using microwave spectroscopy to map out 
the spatial structure of the bosonic states hosted in the 
optical lattice. The proposed techniques can identify 
Mott-insulator states with any site occupation number 
and resolve the spatial distribution of sites with differ- 
ent occupation number. We view this method as inter- 
mediate between less general techniques and high- 
resolution direct imaging ^] which will be necessary for 
long-term development of neutral atom quantum com- 
puting. 

This paper focuses on experimental parameters simi- 
lar to those first used to observe the superfiuid to Mott- 
insulator phase transition in an atomic gas 15]. We 
consider using ^^Rb atoms confined in an optical lat- 
tice formed from three pairs of intersecting laser beams 
with wavelength A — 850 nm. The resulting lattice has 
cubic three-dimensional symmetry with lattice spacing 
a = A/2. The effect of an external potential that changes 
over many lattice sites, formed from magnetic or optical 
fields, is included. The lattice depth is Vb = 30£^r, where 
the recoil energy Ey{, is approximately h x 3.7 kHz {h is 
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Plank's constant). The resulting tunneling strength be- 
tween sites is of order w ^ h x 10 Hz. The on-site inter- 
action U between particles can be changed by changing 
Vq or using a Feshbach resonance. For Vo = SOE'r and 
at zero magnetic bias field, U h x 2.5 kHz. 

The outline of this paper is as follows. In Sec. |nl we 
present a scheme for realizing Mott-insulator configura- 
tions with specified occupation numbers. In Sec. IIIII we 
address the stability of the Mott-insulator configuration, 
in particular, against superfluid order and defect forma- 
tion. In Sec. IIVI we describe experiments that exploit 
microwave transitions between hyperfine states to probe 
the Mott-insulator structure. We conclude with a sum- 
mary of our discussions in SecM 



II. CO-EXISTENCE OF MOTT-INSULATOR 
PHASES 

In this section, we describe an optical lattice in a ra- 
dially symmetric geometry that permits the co-existence 
of concentric Mott-insulating regions with differing oc- 
cupation numbers. We derive the constraints on lattice 
parameters necessary for this co-existence. In contrast 
to previous work 0], we employ counting arguments to 
describe the spatial structure of the Mott-insulator state, 
which allows us to consider a general spatial geome- 
try. Our starting point is the Bose-Hubbard Hamiltonian 
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(1) 



where /i is the chemical potential, V{Yi) is the value of 
the external confining potential at site i of the optical 
lattice (at the point in space), b\ and bi are the boson 
creation and destruction operators at site i, fii = blbi 
is the number operator for bosons on site i, U is the 
on-site interaction energy between two bosons, and W 
represents the tunneling that tends to delocalize bosons, 
with w being the tunneling matrix element for bosons 
between nearest neighbor sites on the lattice. 

In this section, we concentrate on the limit w/U ^ 
where tunneling is negligible and H Hq. (For small 
tunneling w/U, W can be treated as a perturbation; we 
return to this in a later section.) The sites are decoupled 
in this limit and the Hamiltonian, Hq, is diagonalized by 
the set of states with specific occupation number on each 
site \n)i = {bl)^'\0)i/^/ni^.. For a given \n)i, the energy 
is En — ^fln + ^{ri^ — n), where /i is an "effective" 
local chemical potential fl = ^ — V(r). The occupation 
number n that minimizes the energy is determined by 
dEn/dn = = -fl- !l{2n - 1), which gives n = § + 
Since n can take on only integer values, the minimization 
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FIG. 1: The mean-field phase diagram for the Bose-Hubbard 
Hamiltonian at first perturbative order in w/U. Here, jl = 
fi — V{r), as described in the text, and Z is the coordina- 
tion number of the lattice {Z = 6 for the three-dimensional 
cubic lattice under consideration here). The shaded regions 
are Mott insulating with the indicated number of bosons per 
site, n. The white regions are superfluid. Notice that the 
superfluid regions at integer values of jl extend all the way 
to zero tunneling (w/U = 0). While the exact location of 
the boundaries between superfluid and Mott regions at finite 
tunneling is different at higher orders in perturbation theory, 
the qualitative features of the phase diagram persist [l^ . 
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The state with n bosons, which we refer to as the n 
Mott-insulator state, is an incompressible Mott-insulator 
phase. 

With finite but small tunneling [w/U ^ 0), the regions 
where jl/U is nearly an integer become superfluid, lead- 
ing to the well-known phase diagram shown in Fig. ^ 
At finite temperature but negligible tunneling, all Mott- 
insulator phases are replaced by normal fluid phases of 
bosons which display thermally activated hopping 16j. 
At suitably low temperatures kT ^ AE', where /S.E is 
the excitation energy for adding or removing a particle, 
the system is most likely to be found in its ground state, 
which is the situation we consider in this section; we dis- 
cuss the effects of finite temperature and tunneling in 
subsequent sections. 

For a spatially varying external potential V{y) , and 
therefore spatially varying effective chemical potential /i, 
the system can have co-existing regions of different Mott- 
insulator states P, ^3 ■ If the external potential is spher- 
ically symmetric and of the form V{r) = ar^, then jl is 
the largest at r = and decreases for increasing r. Note 
that the arguments in this section can be extended for 
more general potentials lacking spherical symmetry. In 
the w/U = limit, the site occupation number m at the 
center is determined by the condition m—l< [i/lJ < to, 
where the chemical potential /i is determined by the to- 
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FIG. 2: Illustration of the bounding values of r for an 
spherically-symmetric confining potential centered at r = 
which supports shells of Mott-insulator phases. 



tal number of particles N . The number of bosons per 
site changes by one at each radius where [/x — V{y)]/U 
is an integer, leading to a shell structure. The radii 
for the boundaries, as illustrated in Fig[21 are given by 
Rn — {[^J' — nU]/a)^^^ (n = 0, 1, . . . , m— 1), where the n- 
th Mott-insulator state lies between i?„ and Rn-i- This 
shell structure is established by the competition between 
interaction and potential energy. The edge of the occu- 
pied lattice is defined by the outermost boundary of the 
Mott-insulator states where the occupation drops to zero 



All other radii can be expressed in terms of Rq: 
Rn^ Ro[l-n/{fi/U)]^/^. 



(3) 



(4) 



For a given interaction energy U and total number of 
bosons TV, co-existence of a total of m Mott shells will 
occur only for a specific range of the parameter a. We 
find that range as follows: the chemical potential ^ is 
determined via the global constraint 



(5) 



If the spacing between sites in the optical lattice, a, is 
small compared with i?o, we can approximate each Mott 
region by a three-dimensional continuous spherical shell 



with inner radius i?„, outer radius Rn-i and uniform 
density n/a^. In this continuum limit. 
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With some rearrangement, we arrive at the condition 
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1 
(9) 



n=0 



Because m — 1 < /i/C/ < to, we can obtain from Eq. 
the upper and lower bounds on x{^/U) necessary to re- 
alize a shell structure with any desired maximum site oc- 
cupancy TO. For example, for ^ = 2 (a harmonic external 
potential). 



3.83 < 
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An 



i^—j < 9.03, TO = 3. (11) 



As an example, we consider the experimental parame- 
ters detailed in the introduction and N = 10^ atoms. To 
support three Mott-insulator shells in a harmonic trap 
V{r) = ar^ the curvature of the trap needs to satisfy 
the constraint h x 8.8 Rz/fim^ < a < h x 15.6 Hz/^m^. 
For a typical value in this range, for instance, fi/U = 2.5 
(and hence a = h x 12.1 Hz/^m^), Eq. ^ and Eq. Q 
show that the size of the entire Mott-insulator structure is 
given by Rq — 22.8 fiTa, and that the radii separating the 
three Mott shells are Ri = 17.6 /im and R2 = 10.2 /xm. 

In this section, we have shown how the spatial struc- 
ture of the Mott-insulator ground state in the absence 
of tunneling depends on a limited number of experimen- 
tal parameters. Based on counting arguments, we have 
derived a series of constraints on these parameters for 
realizing a shell structure with a specific core occupancy. 



III. STABILITY OF THE MOTT-INSULATOR 
PHASES 

The stability of the Mott structure is relevant to ex- 
periments probing the physics described in section^ Of 
particular practical concern is the elimination of num- 
ber fluctuations when using the Mott-insulator state as a 
fiducial state for neutral atom quantum information ex- 
periments. Several proposals suggest creating the Mott- 
insulator state followed by a purification scheme to pre- 
aie a well-defined state with one atom per lattice site 
0. These schemes cannot eliminate sites that are 
initially vacant and suffer inaccuracy when number fluc- 
tuations become large. 
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Deviations from the nested Mott-insulator shell struc- 
ture will be driven by finite tunneling between lattice 
sites and finite temperature. The dynamic evolution of 
the system is a result of the complex interplay between 
the Bose-Hubbard Hamiltonian, dissipation, and finite 
temperature. While a complete study of stability is be- 
yond the scope of this paper, we summarize the effects of 
tunneling and finite temperature and provide estimates 
relevant to a typical experiment. 

We assume that the lattice parameters described in 
the previous section determine the initial conditions for 
the bosons and that the lattice potential is free of im- 
perfections. In the following sections, we first discuss 
global deviations from the Mott-insulator structure com- 
ing from the tendency to form superfluid regions. We 
show that the effects of superfluidity can be made negli- 
gible in the optical lattice. We then discuss the local fluc- 
tuations caused by "particle" and "hole" (p-h) defects in 
the Mott-insulator structure. We analyze the emergence 
and dynamics of p-h defects at finite temperature and 
tunneling. 



A. Superfluidity 

In principle, even the smallest amount of tunneling can 
alter the Mott-insulator phase. At finite tunneling, the 
natural source of instability is the formation of super- 
fluid regions. As shown in the phase diagram of Fig. ^ 
the system is particularly susceptible to forming regions 
of superfluid at the boundaries between Mott states of 
different occupation number. For a specfic set of values 
for the tunneling w, interaction energy U and external 
potential curvature a, we can estimate the size of the su- 
perfluid regions formed between Mott phases. The Mott 
states remain relatively robust against delocalization if 
the expected thickness of the superfluid shells is less than 
a lattice spacing. 

To estimate the size of the superfluid region at T = 0, 
we invoke a mean-field treatment that decouples lattice 
sites even in the presence of tunneling |18l| . To summarize 
the treatment, the effective Hamiltonian takes the form 



shifted to 



Hmf — 



(12) 



Here b = Zw{h) is related to the superfluid order pa- 
rameter and is obtained self-consistently, with Z being 
the co-ordination number (the number of nearest neigh- 
bors per site). Terms involving b and can be 
treated as a perturbation to the Mott-insulator states 
\n). To second order in b, the perturbed eigenstates at 
a single site are 



= |n) + ^-sa+x/n + l|n + 1) + fsOfi V^|n - 1) 



+*|a+|n + 2) 



(*: 



(13) 



where a^^r^ are functions of n, ^ and U . The energies are 



En = En+Xn\^l 



o 



(l*.l'). 



(14) 



where £"„ is the energy of the unperturbed Mott state \n) , 
and Xn = un\ + ^-uin-i) ■ Self-consistently minimiz- 
ing En, - Zw(b)(b^) - {b)^*B - {b^)^B = En,+p\^B\^ + 

C'^l^s]"'^ determines Here, no is the Mott ground 
state occupation number appropriate for jl. The super- 
fluid boundaries are determined by the sign of p(^, w, U); 
(b) has a non-zero expectation value where p{fi, w, U) is 
negative. We have to first non-vanishing order in Zw/U, 
P = Xno{fJ-/U)[l - ZwXnoifJ-fU)] with 



Xno 



n„oitJ-/U) + 1 

Uno{ii/U)-fi fi-U[noip/U)-iy 



no{fi/U) 



(15) 



The superfluid boundary satisfies the condition 1 



Zwxn 



0. For fixed tunneling, this determines the 



boundary values of the effective chemical potential along 
each Mott-insulator lobe in Fig. ^ For Zw/U ^ 1, 
we have for the effective chemical potential at the up- 
per (/l^g) and lower (/i~jj) boundaries of the hq Mott- 
insulator lobe 



U 



noil- Zw/U), 



U 



{no + Zw/U). (16) 



Ko-i and 



Eq. (|16|) can be used to determine the radii 

between which superffuidity will exist in the con- 
tinuum limit. If the distance between the radii where 
superfluidity can exist is smaller than a lattice spacing 



(A? 



' no-l 



< a), then number fluctuations due 



to superfluid tendencies are expected to be neghgible. 

As an example, we again consider the experimental 
situation detailed in the introduction with a harmonic 
external potential and fi/U = 2.5. The thickness of the 
superfluid between the no = 1 and no — 2 Mott-insulator 
states is Ari_2 ~ 0.82y^U/a Zw/U. Using the value 
a ~ hx 12.1 Mz/ jjLTT? from the previous section, we obtain 
Ari_2 ~ 0.07 /im and Ar2-3 ~ 0.24 /im, which is smaller 
than the lattice spacing a w 0.43 /im. Therefore, for this 
specific case, no superfluid regions will exist. Eq. Hlt)|) 
could also be used to determine the experimental param- 
eters leading to the existence of superfluid shells, which 
may be of fundamental interest in and of themselves. 



B. Defects 

While global instabilities coming from a tendency for 
the bosons to condense can be made negligible, local in- 
stabilities can still permeate the Mott-insulator structure 
in the form of defects. In this section, we discuss two 
types of local excitations — isolated defects where a site 
has an extra particle or one less (hole) compared with 
the Mott-insulator ground state, and excitations where a 
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particle has hopped to a neighboring site and leaves be- 
hind a hole. The former, which we call "site-defects" can 
arise from random removal of atoms from the lattice, for 
example by heating from spontaneous emission driven by 
the lattice light or by collisions with room-temperature 
residual gas molecules in the vacuum system. We refer to 
the latter as "particle-hole" (p-h) excitations and these 
are natural perturbations to the Mott-insulator ground 
state. Site-defects can also arise when components of p-h 
excitations (due to finite temperature or imperfect adi- 
abatic transfer) that were present during initial loading 
of the Bose-Einstein condensate into the lattice become 
widely separated. 

In the following discussion, we first consider the en- 
ergy gap associated with these excitations; the energy 
gap plays a key role in the statics and dynamics of de- 
fects. We then study the static perturbative effects of 
p-h excitations on the Mott-insulator state in the pres- 
ence of finite tunneling. Next, we consider the dynamics 
of defects within the framework of a toy two-site model. 
We comment on defect dynamics over the entire system 
and on the effects of dissipation. Finally, we discuss the 
role of finite temperature. 

Energy Gap In the Mott-insulator ground state at 
zero tunneling, the functional E{n) = —fin + !jn{n — 1) 
is minimized at each site. The energy gap associated with 
single site defects is Ai?+ = Eng+i — Eno = —A + Uni^ 
for the addition of an extra particle and A_B^ = En^-i — 
Eno ~ i^~U{nQ — l) for the removal of a particle (addition 
of a hole) , where hq is the ground-state occupation num- 
ber. Since riQ — 1 < |j < no, AE^ is positive and of order 
U in the bulk of the system. Close to Mott boundaries, 
though, the energy gap becomes arbitrarily small as |j 
approaches either uq or ng — 1. Because of the underlying 
lattice and inhomogeneous confining potential, however, 
the energy cost for a site defect rarely vanishes: there 
is always a change in potential energy AT^ over a lattice 
spacing and it is unlikely for a boundary to exactly coin- 
cide with a lattice site. For the experimental parameters 
detailed in the introduction, Vi = {h x 12.1 IIz//im^)r?, 
we have A^ = hx 181 Hz (h x 104 Hz) at the boundary 
Ri ^ 17.6 fim (i?2 = 10.2 /im). 

The energy cost for a p-h excitation in the bulk is 



E^°P{no) = AE+ + AE-, = AV, + U, 



(17) 



where i is a site index, AVi = Vi — Vi+i, and Vi is the on- 
site potential energy from the confining potential. For 
a particle hopping across a boundary between a Mott- 
insulator no phase on site i and no — 1 phase on site i + 
the energy cost is 

£;f°P(noi) ^ AE-{no) + AE+,,ino-l) = AV, 

E^°^ino T) ^ AE+ino) + AE-^,{no 1) 

= 2U-AV, (18) 

The symbol | denotes hopping from an np to an no — 1 
Mott-insulator state and t denotes the reverse process. 



The most favorable p-h excitation therefore involves a 
particle hopping across a Mott-insulator boundary from 
higher to lower occupation. 
Finite Tunneling - Statics: 

At small tunneling {w <ti U), corrections to the w — 
Mott-insulator ground state jV'o) can be written in terms 
of p-h excitations provided that w <C E^°^ . As found in 
the previous discussion, since the smallest E^°^ is across 
Mott boundaries with E^°^{no I)— AV, the change in 
potential energy must be much larger than the tunneling 
energy (AV^ ^ w), which is achieved by the parameters 
under consideration, namely, w ~ h x 10 Hz and AV = 
h X 181 Hz {h X 104 Hz) across the boundary at i?i — 
17.6 /im (i?2 = 10.2 //m). This inequality is related to the 
condition found in Sec. IIII Al for the putative superfluid 
regions to occupy less than a lattice site. 

To first order in w/E^°p, the perturbative groundstate 
lip) of Eqn.(^ is given by 



l^o) 



E 

cxc 



exc) 



{e^c\W\^o) 
Ai?^™ 



(19) 



where c is a normalization constant. The excitations and 
corresponding activation energies are |exc) and AE'^^'^, 
respectively, and the sum is over all p-h excitations. 
Comparing |-0) and IV'o), we can estimate the deviation 
of the many-body state or the single site occupancy from 
the zero-tunneling Mott-insulator ground state. 

As an example, consider the case of a finite-sized uni- 
form lattice with no external potential and occupancy 
no at each site in the absence of tunneling. To lowest or- 
der, tunneling perturbs the Mott-insulator ground state 
IV'o) = |«o,«o, --na) to 



lip) = jno, ■ 

w^/no{no + l) 
U 



,no) 



(20) 



|no, • • • , (no±l)i, (noTl)i • • • no) 



Deviations of this many-body state from IV^i) can be 
quantified by the zero-temperature fidelity [l3|, defined 



as 



(21) 



1 — /o is essentially the probability that the system is not 
in the zero-temperature Mott-insulator ground state. For 
the wavefunction in Ean. (|20|l . the fidelity takes the form 



/o 



w'^no{no + l)Z 



C/2 



(22) 



where Ng is the total number of occupied sites. A 
high fidelity requires the possibly severe requirement 
\/NgZnQw/U 1. For experiments involving local 
probes, however, the more relevant quantity is the de- 



fect probability, P^^ , at a given site i. The probability 
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for a site to be in the state |no ± 1) can be obtained from 
Eq. 



p: 



dcf 



Z w^no(no-l-l) 



C/2 



C/2 



.(23) 



Low defect probabihty locaUy, then requires 
\/ NsZwtiq/U <C 1, which is satisfied for the exper- 
imental parameters considered in the introduction. 

Modifying these arguments for the Mott-insulator 
shells in an inhomogeneous potential is straightforward. 
The number of p-h excitations varies spatially with the 
energy gap for excitations. The fidelity will therefore 
contain terms of order ISV from sites at boundaries and 
contributions of order U from sites in the Mott-insulator 
bulk. The condition for low defect probability for sites in 
the bulk is the same as for the uniform case. The prob- 
ability of a defect, P^^ , for a boundary site becomes of 
order wno/ IS.V <^ 1, where no is the average number of 
bosons per site near the boundary. 

Finite Tunneling - Dynamics: To analyze the evo- 
lution of a configuration of defects in time and the emer- 
gence of p-h excitations driven by quantum dynamics, we 
use a two-site version of the Bose-Hubbard Hamiltonian: 



H = -w{b\b2 + blbi) + j ^ n,{n,~l) + AVn2. (24) 

1=1,2 

Here, i = 1,2 are site indices and AV is the change in 
on-site potential between sites 1 and 2. 

As the simplest case, we consider the restricted Hilbert 
space of just one particle present, i.e., states |1,0) and 
|0,1) in the |ni,n2) basis. This situation could corre- 
spond site-defects composed of either a hole in the n = 1 
bulk or a particle in the n = bulk, or to the boundary 
between the n — 1 and n — states. For AT^ = 0, cor- 
responding to site-defects at fixed radius in the external 
potential, the quantum state tunnels back and forth be- 
tween |1,0) and |0, 1) at a rate w/h. For AV ^ and 
the particle initially in state |1, 0), the probability for the 
particle to remain in state |1, 0) is 



^Si(t) = l- 



w/AV 



[1 + {2w/AVy 



-(1-cosait), (25) 



where ai = \/ AT^^ -f Aw^ jh. The probability to remain 
in 1 1, 0) is depleted at most by w/AF for AV ^ w, which 
is true for the experimental conditions described in the 
introduction. 

Interactions come into play when there are at least two 
particles and the Hilbert space is confined to ni+n2 — 2. 
To focus on the effect of interactions, we set AV = 0. 
The state |1, 1) corresponds to a state in the n = 1 in 
the Mott-insulator bulk, while |0, 2) and |2, 0) correspond 
to p-h excitations. Only the states |11) and \ip2+) = 
(|20) + |02))/V2 are coupled via the hopping term. The 

(2) 

probability PiJ,i{t) for particles to remain in the state 
1 11) is identical to -PiHi(0 from Eq. H25|l with U replacing 



AV and double the value of w. The probability for |11) to 
decay into p-h excitations is at most of order w/U, which 
is 4 X 10^'^ for the experimental parameters discussed in 
the introduction. The same probability and rate hold for 
the decay of any p-h excitations already present. 

We expect the most pronounced long-range dynamics 
to come from site-defects tunneling along sites of equipo- 
tential. In fact, since the kinetic energy of a site-defect is 
reduced by tunneling, site-defects will tend to delocalize, 
corresponding to superfluid-like behavior. As the sim- 
plest possibility, we can model a particle in the n — 
Mott bulk or a hole in the n = 1 Mott-insulator bulk 
moving along an equipotential surface by setting the po- 
tential energy AV and the interaction energy U to zero 
in the Bose-Hubbard Hamiltonian. In this purely tight- 
binding limit, the amplitude for a particle starting at 
site labelled by an integer vector r to propagate to a site 
labelled by s is given by 



(s|e-^^*|r) 



{2ny 



-/2i:.An-Q^^^^(_2^^) (27) 



where d is the effective spatial dimension, k is the mo- 
mentum, Ar = (r — s), and Jy(a;) is the Bessel function 
of the first kind. The maximum probability amplitude is 
transferred from r to s on a time-scale t h\r — s|/ {2w), 
suggesting that site-defects can exhibit ballistic motion 
along the equipotential surfaces. 

We have so far neglected the effects of dissipation. One 
source of dissipation is spontaneous emission driven by 
the lattice light which will occur at relatively low (1- 
10 Hz) rates. We expect the major source of dissipation 
at any given site to come from the interactive coupling 
of the site to the entire system. The calculation of such 
a dissipative term for optical lattices has not been per- 
formed to our knowledge. Tunneling of site-defects be- 
tween nearest neighbors in the presence of dissipation 
can be treated phenomenologically by the well-studied 
"spin-boson" Hamiltonian, which describes the physics 
of a particle in a double well system coupled to an en- 
vironment We have yet to analyze the effects of 
dissipation on defects in optical lattices, and at present, 
we refer the reader to the exhaustive spin-boson litera- 
ture t;he analogous situation of defects in sohds, 
dissipation arising from coupling to phonons is believed 
to be "superohmic" , leading to weak damping of the os- 
cillations of the equivalent two- level system (lol [lll | . 

In the absence of dissipation, there is no complete tun- 
neling between ground and excited states. However, we 
would expect dissipation to relax an excited metastable 
state to the ground state over a long time scale Tdis = 
[7(2^)^]"^, where 7 is the dissipative rate and AE is 
the energy of the excited state. Dissipation and decoher- 
ence of the site-defect wavefunction can also result in a 
cross-over to classical behavior. In this case, we expect 
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site-defects to propagate diffusively rather than balhsti- 
cally. 

Finite Temperature : The stabiUty of the Mott- 
insulator shell structure requires temperatures suffi- 
ciently low for thermal excitations to be insignificant 
(fcsT ^ U, AV). The state at finite temperature is cap- 
tured by the statistical density matrix 



\En)(En\ 



(28) 



where ks is Boltzmann's constant, T is the temperature, 
and En and \En) are excitation energies and correspond- 
ing states, respectively. The finite temperature general- 
ization of the fidelity given in Eq. 1)21(1 becomes 



/= (V'o|p|V'o)/Tr(p), 



(29) 



where |-0o) is again the Mott-insulator ground state. At 
low temperature, p-h excitations occur with Boltzmann 
weight and the effect of tunneling can be accounted for 
perturbatively. The finite temperature fidelity in Eq. H29|) 
for deviations from the uniform Mott-insulator ground 
state is given by 



fiT) « /o 



1 + ZNsC '=B^ ' 



(30) 



where /o is the fidelity at zero temperature. Likewise, 
the probability of finding a defect at a site i, given by 
Eq. ((23), becomes 



(31) 



Compared to zero temperature, the defect probability 
at finite temperature becomes enhanced by a Boltzmann 
term appropriate for p-h excitations. The defect proba- 
bility is largest across a boundary: 



1 -I- me" 



+ ZNse "bi 



(32) 



where is the number of pairs of sites across the bound- 
ary, AV is the potential difference, and Pf^{0) is the 
corresponding zero-temperature defect probability. 

We have simulated the equilibrium configuration for 
a system of iV « 3 x 10^ atoms using a Metropolis al- 
gorithm, in contrast to a slave-boson model employed 
in [23 ■ Fig. 13 shows the spatial profile of occupation 
number from this simulation for various temperatures. 
The simulation indicates that the temperature must be 
smaller than U /{lOks) for thermal fluctuations of parti- 
cle number to be substantially suppressed. A plot of the 
probability distribution for different site occupancies is 
shown in Fig. 0] Well-defined spatial regions where most 
sites are occupied by the same number of atoms disap- 
pear for temperatures higher than {//(lOfcs). The data 
in Fig. ^c) also show how thermally activated hopping 
is most likely at boundaries between the Mott-insulator 
states closest to the minimum of the external potential. 




— • 
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— •— 
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■--0-- 
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— •— 
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FIG. 3: Spatial profile of occupation number for six different 
temperatures (see legend). A Metropolis algorithm is used 
to calculate the distribution of roughly 3 x 10* atoms among 
lattice sites with a harmonic external potential characterized 
by a = /i X 113 Hz//xm^ and the experimental parameters 
detailed in the introduction. The average number of atoms 
per site (n) is shown at radii r which are multiples of the 
lattice spacing. The radius is measured from the minimum of 
the external potential. 



The time scale for thermally activated hopping can be 
probed in experiments and will play a role in schemes 
where a zero entropy state is initialized for neutral atom 
quantum computing |l4j . At the crudest level of approxi- 
mation, the tunneling rates due to quantum and thermal 
effects can be added to obtain the net rate, particularly 
if the time scales are widely separated [lj| . At zero tem- 
perature we found that dynamic fluctuations of p-h exci- 
tations and site-defect propagation across potential bar- 
riers are largely suppressed due to interactions. However, 
at finite temperature, these processes can be thermally 
activated. To describe hopping between sites, one can in- 
voke a phenomenological double well model and appeal 
to the exhaustive knowledge of reaction-rate theories for 
hopping between metastable states [T^ . 



IV. MICROWAVE SPECTROSCOPY 

In this section we describe a technique for probing 
the spatial structure and energy spectrum of the Mott- 
insulator state discussed in section |n] The method, an 
extension of work in 8] , is compatible with any lattice 
geometry and is capable of resolving the effects of finite 
temperature. In contrast to the previous sections of this 
paper where only one atomic internal state was consid- 
ered, multiple hyperfine states are employed in our spec- 
troscopic technique. Transitions between hyperfine states 
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r(^im) 



FIG. 4: Probability distribution of occupation number for six different temperatures: (a) [//(I00fc_B), (b) {//(SOfcs), (c) 
[//(lOfcs), (d) [//(Sfcs), (e) U/ks, and (f) 2(//fcfl. The probability P„ to find a site occupied by n atoms is shown for radii 
r which are multiples of the lattice spacing. The data are extracted from the Metropolis algorithm used in Fig. |3 and are 
displayed as a stacked column plot (see legend). 



are driven using a microwave-frequency magnetic field. 
Information on the interaction energy U and the density 
profile of the gas is obtained by manipulating the reso- 
nant transition frequency between hyperfine states. De- 
pendence on site occupancy is achieved by changing the 
interaction energy U using a Feshbach resonance, while a 
spatially inhomogeneous, state-dependent optical poten- 
tial provides spatial resolution. 

A. Interaction with an oscillating magnetic field 

We first address the atomic interaction with an os- 
cillating magnetic field in anticipation of explaining the 
detailed experimental procedure in subsequent sections. 
In the presence of an oscillating magnetic field, the exact 
Hamiltonian at one lattice site in the absence of tunneling 
is 

{Hab}+Hl = |^nbK-l) + 

^^flaifla - 1) + Uabflaflb + huJabnb~^ 

-J^^W.B (33) 

i 



where fia and hh are the number operators for atoms in 
hyperfine states \a) and \b), Uaa, Ubb, and Uab are the 
interaction energy for two atoms in states \a) or \b) or 
one atom in each state, the sum runs over the atoms, 
and /i is the atomic magnetic moment operator. The en- 
ergy difference hcoab between states \a) and |6) takes into 
account the effect of any external optical and magnetic 
potentials, including a static magnetic bias field that de- 
fines the z direction. We also assume that the lattice 
potential is identical for both hyperfine states; specifi- 
cally we consider a far off-resonance lattice. The applied 
microwave-frequency magnetic field B = Bi cos [w(t)t] ex 
may have a time-dependent frequency uj. The interaction 
with the magnetic field (the last term in Eq. (|33|l ) can 
be rewritten as Hj = ;in[e*"(*)* + e^^'^^*)*] 
where Ml = —^iabBi/A, fiab is the magnetic moment ma- 
trix element between the \a) and |6) states, and CTj!'' are 
the Pauli raising and lowering operators for atom i. 

Because the states \a) and \b) form an effective spin- 
1 /2 system, the eigenstates of Hab are eigenstates of total 
angular momentum J = J27=i "^'^^ where 5^*^ is the ef- 
fective spin operator for atom i. The interaction with 
the oscillating magnetic field, which can be rewritten 
as Hi = f][e*"(*)* + e'^'^^*)*] ( J+ -I- J_), causes a rota- 



9 



tion of the total spin vector. Since we assume that all 
of the atoms start in one hyperfine state, Hj couples 
states with different rrij, but fixed total angular mo- 
mentum quantum number j = n/2 [n — Ua + Uf, and 
J^\j,mj) = + mj)). All states with j — n/2 

are properly symmetrized with respect to two-particle 
exchange for bosons. 

We take advantage of the Schwinger representation 
and write the eigenstates \j,mj), which are states of 
n spin-1/2 particles, as jn^, In this representa- 
tion, J+\na,nb) = \/na{nb + l)\na - l,nb + 1) and 

J-\na,nb) = y/nb{na + l)\na + l,nb- 1). 

The interaction picture is convenient for calculating 
the effect of the oscillating field. Taking into account 
that the interaction Hamiltonian can only couple states 
where Ua changes by one, that -I- nb is constant, and 
by making the rotating-wave approximation, we obtain 

J_g*['^(t)-(S„„, „,-i)/?i]t I (-34^ 

in the interaction picture. In Eq. 1341 the energy 
En^.nbWa.rib) = HabWa^nb) IS determined by the eigen- 
values of Hab- The quantum state of the atoms on a lat- 
tice site can be written in terms of the time-independent 
eigenvectors of Hab as 



Ua ,nb—n — n. 



\na,nb^n-na) (35) 



(a) 



|0,1> - 
|1,0>- 



(b) 
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FIG. 5: The energy spectrum of Hab for lattice sites with 
n = 1, 2, 3 (a) and the Rabi frequencies for transitions between 
states that are coupled by the interaction Hamiltonian (b) . In 
(a), the energy difference between states and the shift of the 
ground state are shown. 



where n = Ua + rib. In general, the Schrodinger equation 
gives rise to coupled equations for the Cn^.n^ that can be 
solved numerically. 

To illustrate the salient features of this interaction 
Hamiltonian, we will ignore off-resonant excitation and 



discuss the case Lo{t) — En^-i^n^+i — En 



^Sit). While 



for ^^Rb the interaction energy matrix elements are nor- 
mally hyperfine-state-independent to within a few per- 
cent, we consider using a Feshbach resonance so that we 
can have Uab — Uaa and Ubb = PUaa- We discuss only 
the states n = 0,1,2,3 even though the proposed exper- 
imental techniques will work for any site occupancy. For 
typical experimental conditions, three-body recombina- 
tion rates will limit the lifetime of the n > 2 state to a 
few ms. However, we include analysis of that state here 
in light of a recent proposal for reducing three-body 
recombination rates by a factor of 100. 

The energy spectrum of the eigenstates of Hab for 
states with different n is shown in Fig. El^a). Solving the 
Schrodinger equation for the atomic wavefunction gives 
coupled equations of the form 



im (36) 

C'nj,+l,ri5 — 1 — i^na,n,,Cn^,n,,G ^ ^ (37) 



where f2„ 



''na,nb = a/ {na + l)nbfl. In the absence of off- 
resonant excitation and for 6 independent of time, the 



problem reduces to oscillations between two states with 
Rabi frequency fln^,nh (shown in Fig. |Sl[b)). There are 
two important features of the coupled problem that are 
apparent in Fig. [S] The first is that there are no degen- 
erate energy splittings in any manifold for constant n. 
However, the transition frequencies are degenerate when 
the transition involves states with the same Ub between 
manifolds of different n. Also, the Rabi rates for transi- 
tions |0, n) <-> |1, n — 1) and |n, 0) ^ |n — 1, 1) are unique 
to each manifold of constant n, regardless of the value 
of /3 (i.e., |1, 1) ^ |Q, 2) and |2, 1) ^ |1, 2)). In the next 
section, we will describe two experiments that will take 
advantage of these features. 



B. Experimental Proposals 

The two experimental techniques that we describe in 
this section are a simple method for determining the ra- 
tios of sites with different n, and a more developed tech- 
nique for mapping out the spatial structure of the insu- 
lator state. We consider experiments using *^Rb atoms 
and the experimental parameters outlined in the intro- 
duction. With the laser intensity required to reach 30 Ej^ 
lattice depths, spontaneous emission rates driven by the 
lattice light will be negligible on experimental timescales. 
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The limiting timescale for experiments will therefore be 
collisions with residual gas molecules in the vacuum sys- 
tem or spontaneous scattering from an additional optical 
potential. 



1. Rabi oscillations 

The first experiment that we propose uses Rabi oscil- 
lations as a probe of the total population in the lattice 
of sites with any n. A similar technique has been used in 
ion-trap experiments to measure pop ulations in different 
harmonic oscillator levels [2^ l23Ll24j . The total number 
of sites with n atoms is identified by the Rabi oscillation 
amplitude at the Rabi rate f2n,o for the \n, 0) |n — 1, 1) 
transition. Because the resonant frequency and Rabi os- 
cillation rate for this transition are independent of inter- 
action energy, this method is insensitive to the motional 
state of the atoms. 

The gas is first prepared in the Mott-insulator state in 
the lattice; harmonic confinement is provided by an inho- 
mogeneous magnetic field. The atoms in the lattice are 
initialized in the state \a) (through optical pumping, for 
example). A magnetic field at frequency cuab is applied 
to couple the \a) and \b) states, and the population in |6) 
is measured for different lengths of time of the applied 
coupling (via resonance fluorescence, for example). The 
data are fitted to a sum of oscillating functions with fre- 
quencies that differ by a factor ^/k where k = 1,2,3.... 
The total population in the lattice in states of different 
n are then directly connected to the amplitudes of each 
oscillating term in the fit. 

One complication with this scheme is that the inho- 
mogeneous magnetic field will introduce a spread in Uab 
across the lattice, which will cause dephasing of the Rabi 
oscillations at different sites. To suppress shifts in tOat 
across the lattice, the energy difference between states 
I a) and \b) should be insensitive to magnetic field. This 
can be accomplished using a first-order magnetic field 
insensitive transition, for example between the states 
IF = 2,mF = 1) and \F = l^mp = -1) near 3.24 G 
[23 ■ For 10^ atoms, to realize sites with n — 3 will re- 
quire at least a 7 mG spread across the gas using either 
of these states. With a 3.24 G magnetic bias field, there 
will be less than a 27r x 0.1 Hz shift in ujab across the 
occupied sites. 

2. Rapid Adiabatic Passage 

In this section, we explain an experimental technique 
to probe the Mott-insulator spatial structure directly, 
where the occupation of sites is identified using the inter- 
action energy. Rapid adiabatic passage is used to transfer 
atoms in state \a) to state \b). The interaction energy is 
manipulated by using a Feshbach resonance such that 
the resonant microwave frequency shifts for transitions 
between \a) and \b) in a manifold of states with constant 



n. The maximum frequency shift dependence on n per- 
mits unique identification of sites with different n. 

To change the interactions between atoms in state \b), 
we choose \b) = \F = l,mp = 1) and apply a uniform 
magnetic field with magnitude close to the *^Rb Feshbach 
resonance at 1000.7 G jl^ll^l. The magnitude of the field 
is adjusted so that /3 = 5 {Ubb = 5Uaa)- Because this par- 
ticular Feshbach is exceptionally narrow, applying a uni- 
form magnetic field is necessary to avoid spatially- varying 
interaction energy shifts. Therefore, rather than obtain- 
ing spatial information by employing an inhomogeneous 
magnetic field, a spatially inhomogeneous near-resonant 
optical potential is used to shift tOab- The near-resonant 
laser fields are engineered to provide a harmonic exter- 
nal potential; this may be accomplished using diffractive 
optics [2^ I23I to produce a focused laser beam with a 
parabolic intensity profile in three dimensions. To maxi- 
mize the shift in LUab while minimizing spontaneous emis- 
sion rates, we choose |a) = |i^ = 2,mp = 2) and tune 
the laser wavelength and polarization so that the op- 
tical potential for the \b) state vanishes. For different 
atoms, such as *^Rb, with broader Feshbach resonances, 
the complication of employing an optical potential may 
not be necessary. 

The overall effect of the applied, parabolic optical 
potential is to shift LUab — ^^o — 7f ^ parabolically in 
space. The curvature 7 is related to the curvature of 
the intensity profile of the near-resonant optical field, 
luq ~ 9.126 GHz is determined by the atomic hyper- 
fine structure and the magnetic bias field, and r is the 
distance from the minimum of the parabolic optical po- 
tential. The spectrum for transitions between states \a) 
and 1 5) for sites with n = 1,2,3 is shown in Fig. |S1 For 
our experimental procedure, the curvature 7 is adjusted 
so that the minimum interaction energy shift is at least 
equal to the shift from the parabolic optical potential 
for the n = 2 state at the boundary with the n — 1 
state. While this may be incompatible with supporting 
the desired Mott-insulator structure, the applied state- 
dependent potential can be switched quickly before the 
spectroscopy described next is performed. 

The experiment is performed by first preparing the 
atoms in the \a) state. Storing the atoms for long times in 
the 1 6) state while the magnetic field is near the Feshbach 
resonance is avoided to minimize effects due to inelastic 
collisions. A microwave frequency magnetic field is then 
swept from lui to ujh and the number of atoms transferred 
to \b) is measured (using resonance fiuorescence, for ex- 
ample) as the frequency is changed. 

Atoms are transferred to state |6) in stages. As the 
frequency is swept toward ujh, one atom at all occupied 
sites is transferred to \b). This process is interrupted at 
luq by a gap where no atoms are transferred to \b) as the 
frequency is changed. The size of that gap in frequency 
UaaiP — l)/ft — iR\ is related to the interaction energy, 
the size of the region where there are at least two atoms 
per site, and the curvature 7 of the optical potential. 
After this gap in atom transfer, one additional atom per 
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6: The spectrum for transitions between \a) and \b) for 
n = 1,2,3. The transition frequency between states with 
different Ua and rib is shown by the thick hue. The dashed 
hues indicate the boundaries between sites with different n in 
the lattice, which occur at radii _Ro, Ri, and R2- Sites with 
n = 1,2,3 are indicated with green, blue, and red coloring, 
respectively. The interaction energy shifts are indicated, and 
the figure is drawn assuming f3 > 0. 



site is transferred to state \b) for all sites that have at 
least two atoms per site. Another gap is encountered at 
ujQ + UaaiP ~ 1)/^ where no atoms are transferred that is 
2Uaa{P — 1)/^ — 7-^2 wide in frequency. Finally, a third 
atom is transferred to \b) in each site where n = 3 as the 
frequency is swept to uj^. No atoms are transferred to 
|6) as the frequency is increased past ujq + 3Uaa{f3 — 1)/^- 
The interaction energy and the distribution of sites in 
space width different occupancy n can be inferred from 
the number of atoms transferred to \b) as the frequency 
is changed. 

This sequence of transfer between hyperfine states sep- 
arated by gaps is generic to any Mott-insulator structure. 
Transfer begins at w^/ = ~ 7^0' which corresponds 
to the boundary of the occupied sites at Rq. The width 
of each gap where atoms are not transferred to \b) as 
the frequency is changed is kUaaiP ~ T-^L where 

k = 1,2,3.... The highest frequency where atoms is trans- 
ferred is determined by the occupancy m at the core and 

is LOo + mUaaiP - 1)/^- 

The results of a simulation of this experimental tech- 
nique is shown in Fig.[7|for a gas consisting of roughly 3 x 
10'' atoms at six different temperatures. The Metropo- 
lis algorithm is used to prepare the initial distribution 
of atoms among lattice sites at each finite tempera- 
ture. An additional harmonic potential characterized by 




(CO - (Oq) /27I (kHz) 



FIG. 7; Simulation of rapid adiabatic passage experiment for 
roughly 3 x 10* atoms at six different temperatures (see the 
legend, where kb is Boltzmann's constant). The number of 
atoms A^i, in state |fe) is shown as the frequency Urf of an 
oscillating magnetic field is swept. 



a — h X 113 Kz/fim^ is applied to generate a Mott- 
insulator state with a core of rt = 3 sites at zero tem- 
perature. This same curvature 7 = 113 x 27r Hz//i77i^ is 
used for the state-dependent optical potential while the 
oscillating magnetic field is applied. The Schrodinger 
equation is numerically solved on each lattice site using 
= 100 X 2n Hz and = 6.4 x IQ-* x 27r • t Hz. 

We would like to have spectral discrimination between 
sites at different radii, which limits the Rabi frequency f2 
since it sets the effective bandwidth for transfer between 
states. Therefore we chose il so that it is comparable to 
the difference in ujab between sites spaced radially by a 
near the minimum of the optical potential. The sweep 
rate for rt is chosen so that Landau-Zener transitions are 
made with high probability. The data in Fig. [3 would be 
generated using a frequency sweep lasting ~ 1.26 seconds, 
which should be comparable with spontaneous scattering 
times and much faster than loss caused by collisions with 
residual background gas molecules. 

From the simulation results shown in Fig. it is clear 
that the effects of a state-dependent interaction persist to 
relatively high temperature compared with the interac- 
tion energy. The smallest frequency at which atoms are 
transferred to \b) decreases with increasing temperature 
as the size of the gas in the harmonic potential grows. 
In principle, the temperature of the gas can be inferred 
from a fit of data obtained using this technique to the 
simulation results. 
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V. CONCLUDING REMARKS 

The possibility of exploring novel quantum phases of 
matter using ultra-cold atoms trapped in optical lattices 
has piqued the interest of many researchers. In this pa- 
per, we have detailed the necessary conditions for creat- 
ing and probing a shell structure of Mott-insulator phases 
in an optical lattice. Instability of the Mott-insulator 
shell structure arises principally from two effects: finite 
tunneling and finite temperature. Finite tunneling of the 
bosons between neighboring sites leads to superfluid re- 
gions between the shells and particle-hole defects in the 
Mott-insulator regions. We have shown that for reason- 
able experimental parameters, these effects can be made 
negligible. The relevant conditions are both that the 
tunneling parameter must be small compared with the 
on-site repulsion {w/U <C 1) and small compared with 
the typical change in the confining potential V between 



neighboring sites {w/ IS.V <C 1) across the Mott-insulator 
boundary. At finite temperature, while the Mott state 
formally ceases to exist, we have found a reasonable sup- 
pression of defects for temperatures lower than U/{lQkh) 
from Monte Carlo simulations of 3 x 10'* atoms. The ex- 
isting literature on quantum tunneling at finite temper- 
ature provides a starting point for a more detailed ex- 
ploration of these issues in this rich system, for instance, 
the effects of the coupling between a single site defect 
and the rest of the system. In this paper we have also 
proposed an experimental probe that has the capacity to 
spatially resolve the Mott-insulator shell structure. 
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